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Abstract 

After a short review of baryogenesis mechanisms, we focus on the charge transport mech- 
anism at the electroweak scale, effective at strong electroweak phase transitions. Starting 
from the one-loop Schwinger- Dyson equations for fermions coupled to bosons, we present 
a derivation of the relevant kinetic equations in the on-shell and gradient approximations, 
relevant for the thick wall baryogenesis regime. We then discuss the CP- violating source 
from the semiclassical force in the flow term, and compare it with the source arising in the 
collision term of the kinetic equation. Finally, we summarize the results concerning the 
chargino mediated baryogenesis in the Minimal Supersymmetric Standard Model. 
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1 Introduction 



The necessary requirements on dynamical baryogenesis at an epoch of the early Universe are 
provided by the Sakharov conditions: (I) baryon number (B) violation, (II) charge (C) and 
charge-parity (CP) violation and (III) departure from thermal and kinetic equilibrium. The 
Sakharov conditions may be realised at the electroweak transition [1], provided the transi- 
tion is strongly first order. C and CP violation are realised in the standard model (SM) 
for example through the Cabibbo-Kobayashi-Maskawa (CKM) matrix of quarks. At high 
temperatures B is violated through sphaleron transitions, and may be responsible for the 
observed baryon asymmetry today, which is usually expressed as the baryon-to-entropy ra- 
tio [2], ub/s ~ 71b/7ti 7 = 7.0 ±1.5 x 10~ n . This is consistent both, with the nucleosynthesis, 
as determined by the observed D/H ratio, and with recent cosmic microwave background 
observations [3]. 

Over the years a large number of scenarios have been proposed to explain the observed 
matter-antimatter asymmetry of the Universe. Broadly speaking, they can be divided into 
two classes: the models based on the grand unfication of forces, and the models based on the 
electroweak symmetry breaking. In grand unified models [4, 5] one uses the B-L violation 
by the gauge and/or Higgs sectors of the GUT model. The particular realisations include 
baryogenesis at preheating through the inflaton decay into heavy GUT particles [6, 7], whose 
out-of-equilibrium decays can produce a nonzero B-L; leptogenesis [8, 9], which involves 
the physics of heavy majorana neutrinos; the Affleck-Dine mechanism [10], which involves 
the dynamics of (B-L)-violating flat directions of supersymmetric models [11] (which may, 
but need not, be embedded in a grand-unified model), remnants of which may survive as 
B-balls [12]. 

Concerning the electroweak scale baryogenesis, the Standard Model (SM) cannot alone 
be responsible for the observed matter-antimatter asymmetry, primarily because the LEP 
bound on the Higgs mass van > 110 GeV [13] is inconsistent with the requirement that the 
transition be strongly first order, in order for the baryons produced in the symmetric phase 
be not washed-out by the sphaleron transitions in the Higgs ('broken') phase. 

Supersymmetric extensions of the Standard Model, on the other hand, may result in a 
strong first order transition. Indeed, in the Minimal Supersymmetric Standard Model the 
sphaleron bound can be satisfied, provided the stop and the lightest Higgs particles are not 
too heavy, m~ t < 170 GeV and ran < 120 GeV [14]. Non-minimal supersymmetric extensions 
can typically provide stronger phase transitions, and are hence less constrained [15]. 

An efficient mechanism for baryon production at the electroweak phase transition is the 
charge transport mechanism [16], which works as follows: At a first order transition, when 
the Universe supercools, the bubbles of the Higgs phase nucleate and grow. In presence of 
a CP-violating condensate at the bubble interface, as a consequence of collisions of chiral 
fermions with scalar particles in presence of a scalar field condensate, CP-violating currents 
are created and transported into the symmetric phase, where they bias baryon number 
production. The baryons thus produced are transported back into the Higgs phase where 
they are frozen-in. The main unsolved problem of electroweak baryogenesis is the systematic 
computation of the relevant CP-violating currents generated at the bubble interface. Here 
we shall reformulate this problem in terms of calculating CP-violating sources in the kinetic 
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Boltzmann equations for fermions. 

The techniques we report here are relevant for calculation of sources in the limit of thick 
phase boundaries and a weak coupling to the Higgs condensate, which are both generically 
realised in supersymmetric models. In this case one can show that, to linear order in the 
Planck constant H, the quasiparticle picture for fermions survives [17, 18]. In presence of 
a CP-violating condensate there are two types of sources: the semiclassical force in the 
flow term of the kinetic Boltzmann equation, and the collisional sources. The semiclassical 
force was originally introduced for baryogenesis in two-Higgs doublet models [19], and sub- 
sequently adapted to the chargino baryogenesis in the Minimal Supersymmetric Standard 
Model (MSSM) [20]. The semiclassical force corresponds to tree-level interactions with a 
semiclassical background field, and it is universal in that its form is independent on interac- 
tions. The collisional sources, on the other hand, arise when fermions in the loop diagrams 
interact with scalar background fields. These sources arise first from the one-loop diagrams, 
in which fermions interact with a CP-violating scalar background. When viewed in the 
kinetic Boltzmann equation, these processes correspond to tree-level interactions in which 
fermions absorb or emit scalar particles, whilst interacting in a CP-violating manner with 
the scalar background. The precise form of the collisional sources depends on the form of 
the interaction. In the following sections we discuss how one can study the CP-violating 
collisional sources induced by a typical Yukawa interaction term. 



2 Kinetic equations and the quasiparticle picture 

Here we work in the simple model of chiral fermions coupled to a complex scalar field via 
the Yukawa interaction [17, 18] 

£ Yu = -vHl^r - y<i>*4>R*pL, (i) 

which, at a phase transition, may give rise to a complex, spatially varying, mass term 

m(u) = y'& (u) = m R {u) + imj(u) = |m(w)|e ie(u) . (2) 

Here $o( M ) = (^| < ^ ) ( M )|^) ) and 1^) is the physical state. Such a situation is realised, for 
example, by the Higgs field at a first order electroweak phase transition. 

The dynamics of quantum fields can be studied by considering the equations of motion 
arising from the two-particle irreducible (2PI) effective action [21, 22] in the Schwinger- 
Keldysh closed-time-path formalism [23, 24]. This formalism is suitable for studying the 
dynamics of the non-equilibrium fermionic and bosonic two-point functions 

iS aP {u,v) = ip\T c [Mu)^{v)\\n) (3) 
iA(u,v) = (Sl\T c [<l>(u)^(v)]\Sl), (4) 

where the time ordering Tq is along the Schwinger contour shown in figure 1. The C-time 
ordering can be conveniently represented in the Keldysh component formalism. For example, 
for nonequilibrium dynamics of quantum fields the following Wightman propagators are of 
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particular relevance, 
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Figure 1: The closed time contour for the Schwinger-Keldysh nonequilibrium formalism. 

For thick walls, i.e. when the de Broglie wavelength of typical plasma excitations is 
small in comparison to the phase interface thickness L w , it is suitable to work in the Wigner 
representation for the propagators, which corresponds to the Fourier transform with respect 
to the relative coordinate r = u — v , and expand in the gradients of average coordinate 
x = (u + v)/2. This then represents an expansion in powers of £^b/L w . When written in 
this Wigner representation, the kinetic equations for fermions become [25] 



where for simplicity we neglected the contributions from self-energy corrections to the mass 
and the collisional broadening term, which is of the form — e _ *°{S < }{S'/ l }, where Sh = 
(S r + S a )/2 is the hermitean part of the propagator, and S r and S a denote the retarded and 
advanced propagators, respectively. By considering, as an example, the scalar field theory, 
we have been able to show that, to first order in gradients, the flow term can be rewritten 
as [27] 



number, A^ = (A r + A a )/2 and the scalar collision term. This implies that, when working 
to linear order in the width = (i/2)(Il > — and self-energy 11^ = (IF + n a )/2, one 
can include the self-energy and collision term, while the effect of the collisional broadening, 
described by the second term in (7), can be consistently neglected. While we focus here on 
the effects of a pseudoscalar mass, one should keep in mind that one can include the effect of 
the self-energy within the on-shell approximation, provided one appropriately modifies the 
spectral condition. 

When the collision term is approximated at the one-loop (with the resummed propa- 
gators), equation (6) corresponds to the nonequilibrium fermionic Schwinger- Dyson equation 
shown in figure 2. Since the flow term of the scalar equation (also shown in figure 2) does 
not yield CP-violating sources at first order in gradients [17, 27], we do not discuss it here. 

As the bubbles grow large, they tend to become more and more planar. Hence, it suffices 
to consider the limit of a planar phase interface, in which the mass condensate in the wall 
frame becomes a function of one coordinate only, m = m(z). Further, we keep only the 




(6) 




3 




+ 



Figure 2: The one-loop Schwinger-Dyson equations for the out-of-equilibrium fermionic (S) and 
scalar (A) propagators. When projected on-shell and expanded in gradients, these equations reduce 
to the kinetic Boltzmann equations. 

terms that contribute at order h to Eq. (6), which implies that we need to keep second order 
gradients of the mass term 

me -|a- Q k =m+ m 'Q kz _ -m"d 2 k + o(d 3 z ), (8) 
2 8 

where m = m(z), m' = d z m and m" = d 2 m. On the other hand, in the collision term we 
need to consider terms only up to linear order in derivatives 

Ci/j — C^q + Cj/ji + . . 
Ctpo — —~(T; > S < — Y; < S > ^) 

= -^(^^-4 1 >f)(s>5<-S<5>), (9) 

where S < and S > represent the fermionic self-energies, and the derivatives {d^ 2 \ 
d% ) act on the first (second) factor in the parentheses. 

An important observation is that, when G = G(k^, t — xn ■ k», z), the spin 

S z = L-\A)S Z L^) = 7|| (S z -i(v\\ X a) z ) (10) 

is conserved 

[V,S g ]S< = 0, (11) 

where V is the differential operator in Eq. (6), a = 7°7, S z = 7°7 3 7 5 , and 7y = 1/(1 — v 2 ) 1 ^ 2 . 
The spin operator correponds to the boosted spin in the ^-direction (of the interface motion), 
and when written as the Pauli-Lubanski spin operator, 

S PL (k,s) = --^, e = (k 2 y/ 2 , s 2 = -l, s-k = 0, (12) 
e 
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the spin 4-vector corresponds to 



3^ = ^ 



k e 



( k k z \ 
k x k z 
kyk z 

v H J 



(13) 



In the highly relativistic limit we have k^ m k 2 , and the spin vector s oc k, such that the 
spin operator S z approaches the helicity operator, 

H(k) = -— M7 5 
e 

= t/77 5 , h* = ±(] k \.), (14) 

eo y k k J 

as one would expect. (As usually, the helicity operator measures spin in the direction of 

particle's motion, k — k/\k\.) As a consequence, for light particles with momenta of order 
the temperature, k ~ T ^> m, the spin states we consider here can be approximated by the 
helicity states, which are often used in literature for baryogenesis calculations [19, 20]. To 
answer the question to what extent is this fulfilled, requires a detailed quantitative study, 
which is beyond the scope of this talk [27]. 

This discussion implies that, without loss of generality, the fermionic Wigner function 
can be written in the following block-diagonal form 



s< = j:s t 



< 

s 

s=± 



5< = L(A)- 1 ,S < L(A) 
-*T°Sf = \{l + s^)®p a ~gl (15) 

where a 3 and p % {% = 1,2,3) are the Pauli matrices, p° = 1 is the 2x2 unity matrix, and 
L(A) is the following Lorentz boost operator 

L(A) = fc b + ^-7°7-fc|| > (16) 
^2k (k + k ) 

with k = sign(/c )(A;o — kj) 1 ^ 2 , and A corresponds to the Lorentz boost that transforms away 

With the decomposition (15) the trace of the antihermitean part of Eq. (6) can be written 
as the following algebraic constraint equation [18] 

(V-M 2 + f |m|Vk = 0, (17) 

kr 







where gQ Q = ^g^ denotes the particle density on phase space {k^Xv}. Equation (17) has a 
spectral solution 
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2tt 

#00 = ^2-y—n s S{k ^u s ±), (18) 



where u s ± denotes the dispersion relation 



\m\ 2 9' 

u s ± = uj =F s ^ -z , uj = \Jk 2 + \m\ 2 , & = }Jvo-k^ (19) 

and Z s ± = 1 =F s\m\ 2 9' /2lUq. The delta functions in (18) project n s (/c M , t — x\\ - fcy, 2;) on-shell, 
thus yielding the distribution functions f s+ and / s _ for particles and antiparticles with spin 
s, respectively, defined by 

f a+ = n a (u a+ ,k z ,t-x\\'k\\,z) 

— * 

f a - = l-n a (-u 8 -,-k z ,t + xwk\\,z). (20) 

This on-shell projection proves the implicit assumption underlying the semiclassical WKB- 
methods, that the plasma can be described as a collection of single-particle excitations with a 
nontrivial space-dependent dispersion relation. In fact, the decomposition (15), Eq. (17) and 
the subsequent discussion imply that the physical states that correspond to the quasiparticle 
plasma excitations are the eigenstates of the spin operator (10). 

Taking the trace of the Hermitean part of Eq. (6), integrating over the positive and 
negative frequencies and taking account of (18) and (20), one obtains the following on-shell 
kinetic equations [18] 

d t f s ± + v\\ ■ V||/ a± + v s ±d z f s ± + F s ±d k J s± = <Vfc[/ a ±], (21) 

— * 

where f s ± = f s ±(k,z, t — v\\ ■ x\\), C s ±[f s ±] is the collision term obtained by integrating (9) 
over the positive and negative frequencies, respectively, the quasiparticle group velocity 
v s ± = k z /uj s ± is expressed in terms of the kinetic momentum k z and the quasiparticle energy 
u s ± (19), and the semiclassical force 



|m| 2 ± s(\m\ 2 6')' 



In the stationary limit in the wall frame the distribution function simplifies to f s ± = f s ±(k, z). 
When compared with the 1+1 dimensional case [17], the sole, but significant, difference in the 
force (22) is that the CP-violating |m| 2 ^'-term is enhanced by the boost-factor 7y = u /uj , 
ljq = (cUq — fci?) 1 / 2 , which, when integrated over the momenta, leads to an enhancement by 
about a factor two in the CP-violating source from the semiclassical force. 



3 Sources for baryogenesis in the fluid equations 

Fluid transport equations are usually obtained by taking the first two moments of the Boltz- 
mann transport equation (21): integrating (21) over the spatial momenta results in the 
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Figure 3: The flow term sources (25)-(26) characterised by the integrals x\j a {xi) (red solid) and 
xfj'b(xi) (green dashed) as a function of the rescaled mass X{ = \rrii\/T. The sum of the two sources 
(dotted blue) is also shown. 



continuity equation for the vector current, while multiplying by the velocity and integrating 
over the momenta yields the Euler equation. The physical content of these equations can be 
summarized as the particle number and fluid momentum density conservation laws for fluids, 
respectively. This procedure is necessarily approximate simply because the fluid equations 
describe only very roughly the rich momentum dependence described by the distribution 
functions of the Boltzmann equation (22). The fluid equations can be easily reduced to the 
diffusion equation which has so far being used almost exclusively for electroweak baryoge- 
nesis calculations at a first order electroweak phase transition. A useful intermediate step 
in the derivation of the fluid equations is rewriting Eq. (21) for the CP-violating departure 
from equilibrium Sf S i = 5f S i + — 6f S i- as follows 




dk z Sf si + VnSFsiiduft 



+ v w F 0i 5u si (JO-HL) - (dlfu)^ 
where i is the species (flavour) index, f u = \j(e^ w + 1), and 



Cipsii (23) 



Foi 


1 12' 

\mi\ 






(\m t \%)' 
= s 




SF si 


— - 1 Sl+ - 1 SI— 






= ^^si-\- ^-"tpsi — 





(24) 

When integrating (23) over the momenta, the flow term yields two sources in the continuity 
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Figure 4: The semiclassical force baryogenesis mediated by charginos of the MSSM calculated in 
the helicity basis. The figure shows contours for the baryon-to-entropy ratio in the units of 1CP 11 
for two wall velocities v w = 0.01 and v w = 0.03 as a function of the soft susy breaking parameters 
fi and m2- A maximal CP violation in the chargino sector is assumed. The shaded (yellow) regions 
are ruled out by the LEP measurements. The observed baryon asymmetry is in these units 5 — 9. 
(The figure is taken from the latter reference in [20]). 



equation for the vector current. The former comes from the CP- violating spin dependent 
semiclassical force, and has the form 



d 3 k 



= S^w ;V Ja(Xi), 25 

with Xi = \rrii\/T, while the latter comes from the CP-violating shift in the quasiparticle 
energy, and can be written as 

d 3 k 



S", = r, r I 3 F 0i 5u} si 



' 2tt 2 T 



Iml' Jb(xi). (26) 



The total source is simply the sum of the two, S S i = + S^. To get a more quantitative 
understanding of these sources, in figure 3 we plot the integrals J a and in equations (25) 
and (26). A closer inspection of the sources and indicates that the total source S si 
can be also rewritten as the sum of two sources: the source oc |mj| 2 '^', characterized by 
xjJa + x 3 Jb, and the source oc \mi\ 2 9", characterized by x\j a . We note that, in the spin 
state quasiparticle basis the flow term sources appear in the continuity equation for the 
vector current, while in the helicity basis, which is usually used in literature [19, 20], the 
flow term sources appear in the Euler equation. In figure 4 we show recent results of baryo- 
genesis calculations of Ref. [20] based on the CP-violating contribution to the semiclassical 
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force in the chargino sector of the Minimal Supersymmetric Standard Model (MSSM). This 
calculation is based on the quasiparticle helicity states picture, and it suggests that one 
can dynamically obtain baryon production marginally consistent with the observed value, 
n B /s = 7.0 ± 1.5 x 10 -11 , provided m 2 ~ ~ 150 GeV and v w ~ 0.03 (c = 1). 

At this moment, the question of baryogenesis mediated via the charginos of the MSSM is 
not completely resolved. Indeed, more recenly the results, which we summarise in figure 5, 
have been reported [28], where the relevant CP- violating sources were computed in the 
flavour basis. Flavour mixing was, however, not taken account of, which is, on the wall, 
formally of the order h°, and hence cannot be neglected. 

In conclusion we note that, even though we have recently witnessed important progress 
modelling dynamical baryon production in supersymmetric models, some important ques- 
tions remain unresolved which may have quantitative impact on the final results for baryon 
production. 
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Figure 5: Baryon production mediated by charginos of the MSSM from Ref. [28], expressed as 
the multiple of the observed value, t]bbn = tib/s ~ 6 x 10~ n . The figures show contours for the 
baryon-to-photon ratio in the units of 10 -10 as a function of fi and mi for tan f3 = 10, (a) ni2 = // 
and (b) rri2 = 200 GeV. The wall velocity is taken to be v w = 0.01, and the maximum CP violation 
in the chargino sector is assumed. 

We now turn to discussion of the collision term sources in Eqs. (21) and (23). We assume 
that the self-energies S > ' < are approximated by the one-loop expressions (cf. figure 2) 

£<<>(£;,£) = 

d 4 k'd 4 k" 



(27T) 



[(2tt) 4 5(A; -k' + k'^PLS^ik', x)P R A > ' < {k", x) 



+ (2n) 4 5{k — k' — k")P R S < ' > (k',x)P L A < ' > {k",x)], (27) 

where A < and A > denote the bosonic Wigner functions. This expression contains both 
the CP-violating sources and relaxation towards equilibrium. The CP-violating sources can 
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be evaluated by approximating the Wigner functions S >,< and A >5< by the equilibrium 
expressions accurate to first order in derivatives. The results of the investigation are as 
follows. There is no source contributing to the continuity equation, while the source arising 
in the Euler equation is of the form [25] 

n r d A k k ZyJ 2 s\m\ 2 9' 

2 LWY^ = VwV 32^T^ (H ' m ^' (28) 

where the function Xy(|m|,m^) is plotted in figure 6. It is encouraging that the source 
vanishes for small values of the mass parameters, which suggests that the expansion in 
gradients we used here yields the dominant sources. Note that the source is nonvanishing only 
in the kinematically allowed region, > 2\m\. When the masses are large, \m\, ^> T, the 
source is, as expected, Boltzmann-suppressed. It would be of interest to make a comparison 
between the sources in the flow term and those in the collision term, and apply our methods 
to realistic models [27]. 
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Figure 6: The collisional source contributing to the fermionic kinetic equation at one loop for the 
mass ratios m^j\m\ = 2.1,2.5,3,4,10 and 20, respectively. The source peaks for \m\ ~ 0.7T and 
777,0 w 4|m|. 
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